ALGEBRAIC MEAN FIELD THEORY 



Ts. Dankova and G. Rosensteel 
Physics Department, Tulane University, 
New Orleans, Louisiana 70118 

February 9, 2008 

Abstract 

Mean field theory has an unexpected group theoretic mathemat- 
ical foundation. Instead of representation theory which applies to 
most group theoretic quantum models, Hartree-Fock and Hartree- 
Fock-Bogoliubov have been formulated in terms of coadjoint orbits 
for the groups U(n) and 0(2n). The general theory of mean fields is 
formulated for any arbitrary Lie algebra g of fermion operators. The 
moment map provides the correspondence between the Hilbert space 
of microscopic wave functions and the dual space g* of densities. The 
coadjoint orbits of the group in the dual space are phase spaces on 
which time-dependent mean field theory is equivalent to a classical 
Hamiltonian dynamical system. Indeed it forms a finite-dimensional 
Lax system. The SU(3) mean field theory is constructed explicitly in 
the coadjoint orbit framework. 



1 Introduction 

There are two basic types of models widely used in describing an 
atomic nucleus: shell model and mean field theory. These theories are 
the foundations of many-body science. Both of them can be built on 
a Lie algebra of fermion operators with the differences arising from 
the way the algebra determines the model states. In a shell model 
the states of the many-fermion system form an irreducible unitary 
representation of the Lie algebra. The Hamiltonian is an element of 
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the enveloping algebra and its eigenvalues form the energy spectrum. 
The matrix elements of the Lie algebra generators give the transition 
rates between the states. In a mean field theory the model states are 
restricted to an orbit of the corresponding Lie group in Fock space. 
This orbit is a general surface and it is not a vector space. There is 
a well-defined mapping from the orbit to the space of densities. Its 
image is a phase space and time-dependent mean field theory is a 
classical dynamical system on it. 

A familiar illustration of this idea is the algebra u(n) of one body 
operators X = J2 aa ' ^aa' a a a a'^ where (V a , a a i are the fermion creation 
and annihilation operators and a, a' = 1, . . . , n index an orthonormal 
basis for the n-dimensional single particle space. The Hilbert space of 
an A-fermion system is the totally antisymmetric irreducible unitary 
representation of u{n) of dimension n\/A\(n — A)\. The shell model 
states are eigenfunctions of the Hamiltonian 

H = ^ T aa ia) a a a i + — ^ V ai3a> ■p'a ] a a ] l3 apia a i . (1) 

aa' af3a'/3' 

In the case of Hartree-Fock mean field theory, the many-body 
states are approximated by Slater determinants that form an orbit 
of the group U (n) in Fock space. [l[ The group of unitary basis trans- 
formations for the single particle wave function space is U(n). A 
determinant ipHF is typically a poor approximation to the shell model 
eigenstate ip. Indeed these two wave functions are almost orthogonal, 
(ip | ^Phf) ~ 0. But the expectation values of the one-body operators 
are similar whether taken with respect to the shell model ground state 
or its mean field approximation, 

U> | X | V) a (iPhf I X I te), (2) 

for any one-body operator X. 

The fully antisymmetrized wave functions are in 1-1 correspon- 
dence with the idempotent (p 2 = p) density matrices p a i a = (ipHF \ 
($ a a a t | iPhf)- The idempotent densities form a surface in the space 
of all density matrices that is invariant with respect to the coadjoint 
group action, p i— ► gpg~ 1 for g G U (n). The dimension of this coadjoint 
orbit surface is 2A{n — A) for the case of determinants - a more man- 
ageable number than the factorially-growing dimension of the shell 
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model space. The Hartree-Fock Hamiltonian 

HhF = ^2 ^aa' + Vaf3a'/3'Pf3'f3 <^ a a a' (3) 
aa' \ 130' J 

is an element of the algebra u(n) and it depends on the density. 

The expectation of any one-body operator is given in terms of the 
density matrix by (X) = tr{pX). Indeed the density matrix extracts 
from the microscopic wave function only the information about its 
single-particle degrees of freedom. Thus, the Hartree-Fock approxi- 
mation is adapted to provide the optimal description for the single- 
particle degrees of freedom of a many-fermion system. This is achieved 
mathematically by using the Lie algebra u(n) of one-body operators. 
In a similar way Hartree-Fock-Bogoliubov theory has been presented 
as a coadjoint orbit of the group 0{2n) of quasiparticle transforma- 
tions. § 

When other degrees of freedom dominate in an interacting many- 
body system, the method of coadjoint orbits should be generalized 
to the Lie algebra g corresponding to the physically relevant degrees 
of freedom. In this article the algebraic theory of mean fields is for- 
mulated and applied to the Elliott SU (3) algebra. A synopsis of the 
generalization is as follows: The moment map establishes the connec- 
tion between the Hilbert space of microscopic wave functions and the 
state densities. The densities are elements of the dual space g* to the 
Lie algebra. A generalized Hohenberg-Kohn theorem shows that there 
is an energy functional on the space of densities that is inherited from 
the microscopic Hamiltonian. Stationary densities can be obtained by 
minimizing the energy functional. Each coadjoint orbit is an even di- 
mensional phase space that has a nondegenerate symplectic structure. 
The symplectic form enables the mean field Hamiltonian, an element 
of the Lie algebra g, to be derived from the energy functional. Time- 
dependent mean field dynamics on the space of densities is a classical 
dynamical system of Lax type. 

2 Moment Map and Coadjoint Orbits 

Denote the Hilbert space of ^4-fermion wave functions by H. Sup- 
pose a matrix Lie algebra g is represented faithfully by an algebra 
of hermitian fermion operators; the matrix X € g is represented by 
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the hermitian operator X acting on H. Suppose G is the connected 
Lie group corresponding to the algebra g and U is the corresponding 
faithful unitary representation of G on Fock space. The dual space g* 
is, by definition, the vector space of all real-valued linear functionals 
defined on g. The relation between the algebra and its dual space is 
similar to the one between contravariant and covariant vector spaces 
or between ket and bra spaces. In the context of mean field theory the 
dual elements are the state densities and the notation p will be used 
for them throughout this exposition. The bracket notation, (p,X), 
is used for the operation of the linear functional p G g* on the Lie 
algebra element leg. 

The moment map is defined as a mapping M : H — > g* from the 
Hilbert space to the dual of the Lie algebra that provides a unique 
density p = M(ip) for every wave function i|)GH, 

for all X in g. Although the correspondence between wave functions 
and densities is one-to-one for the special case of Hartree-Fock, the 
moment map is many-to-one in the general case. 

The moment mapping is the heart of mean field theory. The den- 
sity retains only part of the entire information about the system that 
the wave function carries, but a very important part - the expecta- 
tions of the observables that close the Lie algebra. Hence the moment 
map provides a significant simplification of the problem of how to de- 
scribe a many-fermion system. It reduces all the degrees of freedom 
incorporated in the wave function to just those relevant to the observ- 
ables of interest. If the algebra g contained all (not just one-body) 
operators, there would be no such reduction, the mapping would be 
one-to-one, and the shell model would be physically and mathemati- 
cally equivalent to the mean field theory. 

Another important feature of the moment map is that it induces an 
action of the Lie group on the dual space of densities from the unitary 
group representation in the Hilbert space. Given any wave function ip 
and its corresponding density M(ip), the group G transforms the wave 
function into TJ(g)ifj while the density is transformed into M(JJ(g)ip). 
But the transformed density is given by 
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_ (jj I u(gr^u(ff) | f, 

WW) 

= (Mty),g- 1 Xg), (5) 

since the matrix Ad g -i(X) = g~ 1 Xg is represented by the hermitian 
operator TJ (g)^ 1 XXJ (g) . By definition, the coadjoint group action 
Ad* on the space of densities is given by {Ad g * p , X) = (p , Ad g -iX). 
Therefore it has been proven that Ad g * (M{ip)) = M(U(g)ip). 

The mean field approximation restricts the allowed densities in the 
theory to a single coadjoint orbit. Fix a density p E g*. The set of all 
densities that can be obtained by acting with the group in all possible 
ways on a given fixed density p is its coadjoint orbit, 

O p = {Ad g * p € g* | g G G} . (6) 

How does one work with coadjoint orbits? They are multidimensional 
surfaces and defining coordinates on them even locally is an extremely 
difficult task. Fortunately, there is a coordinate-free way. Coadjoint 
orbits are diffeomorphic to homogeneous spaces, i.e. they can be pre- 
sented as coset spaces of the group modulo its isotropy subgroup. By 
definition the isotropy subgroup is the set of all group elements that 
fix a point p G g* , 

G p = {geG\Ad* g p = p}. (7) 

Then the coadjoint orbit is diffeomorphic to the homogeneous space 

O p = G/G pi (8) 

and its dimension is less than the dimension of the dual space itself, 
dim O p = dimC — dimG p = dim g* — dimG p . This dimensional 
reduction is one more advantage to working with the coadjoint orbit 
mean field method. In Hartree-Fock the isotropy subgroup is the direct 
product U(A) x U(n — A) and the dimension of the coadjoint orbit of 
idempotent densities is n 2 — (A 2 + (n — A) 2 ) = 2A(n - A). 

3 Mean Field Dynamics 

The ground state density is found by minimizing the energy functional 
defined on the coadjoint orbit O p . This energy functional should be 
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derived from the microscopic Hamiltonian H on H. Because the mo- 
ment map is many-to-one, the energy functional on the dual space 
cannot be set equal to the expectation of the Hamiltonian but instead 
is given by 

E[p]= inf (9) 

The mean field Hamiltonian, an element of the Lie algebra g, is 
determined by both the energy functional and the symplectic structure 
on the coadjoint orbit. The symplectic geometry is inherited from the 
Lie algebra as follows: Every Lie algebra element X defines a tangent 
vector to the coadjoint orbit surface through the density p. This vector 
is the tangent to the curve "fx(t) = Ad*(p) € O p where g = exp(tX). 
Note that the elements of the isotropy subalgebra gp correspond to 
null tangent vectors, and the tangent space is isomorphic to the coset 
vector space g/gp. For any two tangent vectors X, Y £ g/gp to the 
coadjoint orbit through p, the bilinear form, 

u p (X,Y) = (p,[X,Y]), (10) 

is antisymmetric, closed, and nondegenerate. Thus each coadjoint 
orbit is a symplectic manifold. || In particular, each surface O p is 
even dimensional. Moreover the group G acts on each orbit as a 
transitive group of canonical transformations, 

u Ad * {p) {Ad g {X),Ad g {Y)) =u p (X,Y), (11) 

where X, Y are tangent vectors at p and Ad g {X) and Ad g (Y) are 
tangent vectors at Ad*(p). 

The model group theoretic Hamiltonian h may be defined now. 
Choose a coadjoint orbit Op in the dual space and let E denote the 
energy functional for this orbit. If X is a tangent vector at p, denote 
the derivative of the energy in the X direction by 

dE p {X) = jE(Ad* g (p))\ t=0 , where g = exp(iX). (12) 

The model Hamiltonian h € g/gp is determined by 

dE p (X) = uj p {h, X) for all X G g. (13) 

There is a unique solution to this equation for h G g/gp because 
the symplectic form u is nondegenerate on the tangent space. The 
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model Hamiltonian h depends on the point p in the coadjoint orbit, 
as expected. The time evolution of a density on the coadjoint orbit 
surface can be shown to obey the Lax equation j|] 

P=[h(p),p]. (14) 

It can be easily shown that all the constants of motion for a Lax system 
are given by 

In = ~tr{p n ). (15) 
n 

In the special case when the Hamiltonian H is an element of the 
Lie algebra g, the energy functional is just E[p] = (p, H) and the mean 
field Hamiltonian h is equal to H. For Hartree-Fock, the model Hamil- 
tonian h on the coadjoint orbit of idempotent densities is identical to 
the usual mean field Hamiltonian. Otherwise the model Hamiltonian 
h is the optimal approximation to the microscopic H relative to the 
Lie algebra g in the sense that (ip \ [H — h, X] ip) =0 for all X £ g 
where M(ip) = p. 

4 Mean Field Theory of the 577(3) Model 

The angular momentum L and the Elliott quadrupole operator Q close 
under commutation to form the Lie algebra su{3). ||] Its irreducible 
representations are labeled by the pair of nonnegative integers (A , p), 
where A = and n x , n y , n z are the oscillator 

quanta in the x, y, and z directions of the highest weight vector ?Phw- 
The Hamiltonian of the model is 

H su{3) =H + (1/2)kQ.Q (16) 

with being the oscillator Hamiltonian and k < 0. The energies in 
the spectrum are proportional to the eigenvalues of L 2 . 

The Elliott su(3) algebra is a faithful representation of the complex 
antihermitian traceless 3x3 matrices, 

su(3) = {e g m 3 (c) | e f = -e , tr e = 0} . (17) 

The real antisymmetric matrices correspond to the angular momen- 
tum and the pure imaginary symmetric matrices correspond to the 
quadrupole operator. Because su(3) is a semisimple algebra, its Killing 
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form is nondegenerate ( it can be shown that (p, £) = tr(p£) ), and 
therefore its dual space may be identified with the algebra itself 

su(3)* = {p£ M 3 (C) | p f = -p, trp = 0}. (18) 

If the density is expressed as a sum of real and imaginary parts, p = 
A + iB, then (p, Ly) = Aij and (p, Q^) = B{j . 

The coadjoint action of the group SU(3) on the dual space su(3)* 
is Ad* p = gpg^ 1 , and the coadjoint orbit at each point p £ su(3)* is 

O p = {A£ p = gpg- 1 | 5 e 5C/(3)} . (19) 

Each orbit in this case contains a diagonal matrix with pure imaginary 
entries and can be written in terms of (A, p) as follows 



PHW = -5 -A + /i (20) 




Note that this diagonal density corresponds, via the moment map, to 
the highest weight vector ipnw- The coadjoint orbit is diffeomorphic 
to the coset space G/G p , and the isotropy subgroup depends on the 
degeneracies of phw'- 

( SU(3) /SU{3) = {p} , (A , p) = (0 , 0) , dim = 

O p = \ SU(3)/U(2) , (A ^ , fi = 0) or (A = , p ^ 0) , dim = 4 

[SU(3)/U(1) xU{l), (A/0,/j^0), dim = 6 

(21) 

Consider the quadrupole-quadrupole interaction Q.Q. The en- 
ergy functional that corresponds to this interaction is E[p] = trB 2 . 
Using the symplectic structure the mean field approximation to the 
quadrupole-quadrupole interaction is, according to Eq. (|l3l) , h[p] = 
2p t , and the time evolution of the mean field is governed by the Lax 
equation 

p=[2p\p\. (22) 

For a Lax equation, it is easily shown that the trace of any power 
of p is a constant of the motion. For the su(3) highest weight orbit, 
the conserved trace of the square of the density is evaluated at the 
highest weight density to be I2 = 1/2 tr p 2 = (A 2 + Xp + p 2 )/3. At an 
arbitrary orbit point p = A + iB, the quadratic Lax invariant is given 
by/ 2 = l/4L 2 + l/12Q-Q. 
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